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ABSTRACT 


Pressure,  .volume,  and  energy  equation  of  state  data  obtained  using 
shock  wave  techniques  are  presented  for  aluminum  and  Teflon.  Solid 
aluminum  samples  initially  at  room  temperature  or  preheated  to  near 
melting,  and  porous  aluminum  samples  at  room  temperature,  were  studied 
over  a  pressure  range  of  200  to  1200  kfcar.  It  was  found  that  the 
largest  variations  of  volume  and  energy  could  be  achieved  using  porous 
samples.  Values  of  Gruneisen's  ratio  estimated  from  values  of  the 
thermal  pressure  and  thermal  energy  range  from  2.1  to  1.37.  Due  to 
the  sensitivity  of  Gruneisen's  ratio  to  the  Hugoniot  data,  it  is  not 
possible  to  formulate  its  energy  or  volume  dependence  conclusively  at 
the  present  time.  Of  greatest  significance  is  the  fact  that  it  does 
not  vary  widely. 

Solid  and  porous  Teflon  samples  were  studied  over  a  pressure 
range  of  100  to  500  kbar.  Hugoniot  curves  drawn  on  the  basis  of  the 
ten  data  points  obtained  indicate  a  variation  of  Gruneisen’s  ratio  from 
about  0.7  to  2.0. 
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SECT  I  OS  1 


INTRODUCTION 

Calculations  describing  wave  propagation  in  a  continuous  medium  re¬ 
quire  the  knowledge  of  a  constitutive  relation  or  equation  of  state 
characterizing  the  material  properties  of  the  medium.  In  solids  at 
pressures  greatly  in  excess  of  the  shear  yield  stress,  an  equation  of 
state  relating  pressure,  volume,  and  internal  energy  sufficiently  de¬ 
scribes  the  medium.  The  uork  reported  here  was  undertaken  to  study  the 
equation  of  state  of  aluminum  and  Teflon  using  shock-wave  techniques. 

The  largest  portion  of  the  effort  was  spent  investigating  aluminum. 

In  order  to  reach  thermodynamic  states  lying  in  a  region  of  the  pressure- 
volume  plane  rather  than  along  a  single  Hugoniot  curve,  the  initial  state 
of  the  material  must  be  varied.  This  variation  was  achieved  by  using 
solid  and  porous  aluminum  at  room  temperature  and  by  preheating  solid 
aluminum  to  several  hundred  degrees  prior  to  shocking.  It  was  found 
that  the  largest  changes  in  the  final  shock  state  could  be  produced 
using  porous  samples. 

The  following  sections  include  a  general  discussion  of  equations 
of  state  of  simple  systems  and  how  they  may  be  studied  by  shock  waves, 
a  description  of  the  experimental  techniques  used  and  the  performance 
of  the  experiments,  and  an  analysis  of  the  data  and  comparison  with 
other  data  and  equations  of  state. 
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SECT! OX  II 


BASIC  THERMODYNAMICS  OF  SHOCK  WAVES 

1.  EOLATIONS  OF  STATE 

The  term  “equations  of  state”  has  unfortunately  been  used  in  shock 
wave  literature  to  denote  a  variety  of  relationships  which  convey  dif¬ 
ferent  amounts  of  thermodynamic  information  about  the  system  under  con¬ 
sideration.  Quite  frequently  curves  such  as  isentropes.  adiabats,  or 
Hugoniots  are  referred  to  as  equations  of  state.  The  specification  of 
one  of  these  curves  is  not  a  specification  of  the  equation  of  state  of 
the  system  because  an  equation  of  slate  is  a  function  of  at  least  two 
independent  variables.  The  purpose  of  this  section  is  to  summarize  the 
basic  thermodynamics  of  equations  of  state  of  simple,  nonreacting,  single 
ptiase  systems. 

In  a  simple  system  the  two  independent  variables  are  generally  taken 
to  be  any  two  of  the  four  variables  volume,  pressure,  temperature,  and 
entropy,  denoted  by  l'",  P,  T,  and  S,  respectively.  Expressions  for  the 
remaining  two  as  functions  of  the  two  so  chosen  constitute  a  complete 
equation  of  state,  since  they  yield  a  complete  thermodynamic  description 
of  the  system.  Each  of  the  relationships  independently  is  called  an 
incomplete  equation  of  state,  or  simply  equation  of  state.  Given  the 
two  relationships,  T  -  T{S,V),  P  =  P(S,V),  which  can  in  principle  be 
inverted  so  that  any  pair  of  variables  is  independent,  the  internal  energy 
is  given  by 

dE  =  TdS  -  PdV  (1) 

Knowledge  of  the  internal  energy  E(S,Y)  is  then  equivalent  to  the  two 
equations  of  state  T(S,F)  and  P{5,D  since 


T(S,V)  =  j 

If)  ■  e(s,v)  ■  -k 

l) 

\ 

~S/y  \- 

1  's 

The  function  E(S,Y )  is  frequently  called  a  fundamental  equation  or  complete 
equation  of  state  because  it  also  contains  a  complete  thermodynamic 
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then  S(V,T)  and  P(V,T)  form  a  complete 
Eq.  (1),  the  Helmholtz  free  energy  is 


are  chosen 
equation  of 
given  by 


as  independent  variables 
state.  Analogous  to 


dF  =  -SdT  -  PdV  (3) 

and  F(V,T)  is  the  fundamental  equation  equivalent  to  the  two  equations 
of  state  S(VtT)  and  P(VtT)  since 


S(V,T) 


P(V,T) 


(4) 


Similarly  the  enthalpy  function  H(S,P)  and  the  Gibbs  free  energy  G(P,T) 
serve  as  fundamental  equations  when  the  independent  variables  are  the 
pairs  ( S,P )  and  (P,T)  respectively.  Complete  equations  of  state  may  be 
formulated  in  other  ways;  all  formulations  of  a  complete  equation  of 
state  are  equivalent  and  can  be  deduced  one  from  another  through  general 
thermodynamic  relationships. 

Quite  frequently  one  of  the  incomplete  equations  of  state  is  the 
internal  energy  E(V,T)  or  E{P,V).  Recall  that  the  internal  energy  is  a 
complete  equation  of  state  or  fundamental  equation  only  when  it  is  ex¬ 
pressed  as  a  function  of  S  and  V.  If  E(V,T)  is  known,  supplementing  it 
with  either  P{V,T)  or  S{V,T)  constitutes  a  complete  equation  of  state1. 
Similarly  supplementing  E{P,V)  with  either  T(P,V)  or  S(P,K)  constitutes 
a  complete  equation  of  state1. 

In  shock  wave  research  the  incomplete  equation  of  state  E(P,V)  is 
of  particular  interest.  This  equation  of  state  involves  only  the  mechani¬ 
cal  thermodynamic  variables  P,  V,  and  E  and  does  not  yield  values  of  the 
thermal  thermodynamic  variables  S  and  T.  Shock  wave  measurements  yield 
an  E{P,V)  incomplete  equation  of  state  since  the  values  of  pressure, 
volume,  and  energy  are  computed  from  the  Rankine-Hugoniot  jump  conditions 
expressing  the  conservation  of  mass,  momentum,  and  energy,  that  do  not 
involve  the  entropy  and  temperature  explicitly. 


If  one  is  known,  the  other  csn  be  determined  in  principle  from  Eq.  (1). 


•2Vr- 


1!  •  STOfSTSCT- 
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Let  us  then  summarize.  In  a  simple  thermodynamic  system  there  are 
two  independent  variables.  These  two  are  generally  chosen  from  F,  T,  S, 
and  P.  An  equation  expressing  one  of  the  remaining  pair  as  functions  of 
the  independent  variables  is  called  an  incomplete  equation  of  state. 
Expression  of  both  remaining  variables  in  terms  of  the  two  independent 
variables  is  a  complete  equation  of  state,  t.e .,  a  set  of  relationships 
from  which  all  other  thermodynamic  properties  of  the  system  can  in  prin¬ 
ciple  be  computed.  The  internal  energy  may  be  added  to  the  above  quartet 
of  variables,  but  it  is  related  to  them  through  Eq.  (1).  An  isentrope  is 
the  intersection  of  the  incomplete  equation  of  state  surface  P  =  P(S,V ) 
with  a  plane  S  =  constant.  Thus,  prescribing  an  isentrope  P  -  P{SQ,V)  -f(V) 
is  not  prescribing  an  equation  of  state.  Similarly  a  Hugoniot  curve  is 
obtained  by  subjecting  an  E(P,V)  incomplete  equation  of  state  to  the  con¬ 
straint 

E(P,V)  =  £„  +  |  <P„  +  P)(V„  -  V) 

where  £Q  =  E(P0,V0).  Specifying  a  Hugoniot  curve  is  not  specifying  an 
equation  of  state  of  the  system. 

2.  MIE-GRUNEISEN  EQUATION  OF  STATE 

It  was  stated  earlier  that  equation  of  state  information  obtained 
from  shock  wave  data  leads  to  an  incomplete  equation  of  state  of  the 
type  E(P,V)  This  equation  of  state  tells  the  value  of  the  internal 
energy  of  the  system  over  a  region  of  the  pressure-volume  plane  but  does 
not  contain  the  information  necessary  to  compute  the  temperature  or 
entropy  in  the  same  region.  Since  temperature  and  entropy  are  thermal 
thermodynamic  variables,  as  contrasted  to  pressure  and  volume  which  are 
mechanical  thermodynamic  variables,  the  E{P,V)  equation  of  state  might 
be  suitably  called  a  mechanical  equation  of  state.  The  energy,  pressure, 
and  volume  are  obtained  from  shock  data  by  measuring  wave  velocities  and 
applying  the  Rank ine-Hugoni ot  jump  conditions  expressing  conservation  of 
mass,  momentum,  arid  energy. 

In  problems  of  flow  calculations  the  incomplete  E{P,V)  equation  of 
state  may  be  sufficient  since  it  does  permit  calculation  of  the  curves 
of  constant  entropy  in  the  pressure-volume  plane.  The  constant  value  of 
the  entropy  associated  with  each  isentrope  is  unknown,  however.  One 
P-V-E  equation  of  state  which  has  been  extensively  used  in  describing 


m 

0 


i  (.  Mi' 


metals  is  the  Mie-Grune isen  equation  of  state  with  the  GrOneisen  ratio 
r  assumed  to  be  a  function  of  volume  only.  In  general  the  GrOneisen 
parameter,  as  for  all  thermodynamic  state  variables,  is  a  function  of  two 
independent  variables  and  is  defined  by 


/BP\  (BP/BD„  (36'/3lOr 

r(X%Y>  =  VW)v  =  ^  WWTV  =  ^  T(3S/3r)y  (5) 

where  X  and  Y  represent  any  two  of  the  five  variables  E,  P,  S,  V,  and  T 
which  are  chosen  as  independent.  We  shall  define  the  generalised  Mie- 
GrOneisen  equation  of  state  as  the  solution  of  any  one  of  the  partial 
differential  equations  of  Eq.  (5)  when  r(Jf,  K)  is  specified  All  of  the 
definitions  of  Eq.  (5)  are  equivalent  and  may  be  derived  one  from  ^nother 
by  applying  thermodynamic  identities. 

Using  the  last  definition  of  Eq.  (5), 


Along  an  isentrope  S(VfT )  =  constant  and 


Therefore  the  isentropes  in  the  V-T  plane  are  given  by  the  integrals  of 


3T\  _  r~(K,  T) 

*V/s  "  v 


(8) 


The  usual  Mie-GrQneisen  equation  of  state  as  applied  to  solids  is  obtained 
by  writing  the  pressure  and  internal  energy  each  as  a  sum  of  0°K  corapres- 
sional  component  and  a  thermal  component: 


E(V.T)  =  Ec(V)  +  Eth{V,T) 

P(V,T )  P  (K)  +  P.AV.T)  (9) 

C  t  n 

where  f>c(V)  =  ~dE c(  V)/dV,  the  subscript  c  refers  to  0°K.  and  the  thermal 
components  vanish  for  all  V  at  T  -  0.  Also  the  Gruneisen  parameter  P  is 
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taken  to  be  a  function  of  volume  alone,  i.e., 
integration  of  Eq.  (8)  yields 


!"  =  r(l/).  In  this  case 


T  =  WS)9(V) 


where 


fP(V) 

-1—dV 


6{V)  =  80e 


and  <l>(5)  is  function  of  the  entropy  and  hence  constant  along  an  isentrope. 
The  thermodynamic  definition  of  the  temperature  combined  with  Eq.  (10) 


gives 


m,  i')  = 


•  (a 


=  *(S)0(K)  . 


Integration  of  Eq.  (12)  then  yields 


£■(5,10  -  Ec{V) 


=  6 (10  /S«&(S)dS 


since  E(0,  V)  =  Ec(V).  Equation  (13)  is  a  fundamental  equation  since 
it  expresses  the  internal  energy  E  as  a  function  of  S  and  V.  This  equa¬ 
tion  can  be  evaluated  explicitly  as  a  function  of  S  and  V  if  the  functions 
£  (l 0,  $(5),  and  6{V )  or  r(IO  are  known.  Thus,  for  a  system  in  which  T  is 
a  function  of  volume  only,  the  specification  of  Ec(V),  r(K)»  and  4>(S)  con¬ 
stitutes  a  complete  equation  of  state.  The  pressure  is  given  by 


PiS.V ) 


dEc{V)  ddi,V) 


4>(S)dS 


But  ~'dEc(V)/dV  =  Pc(V)  and  therefore 


F(V)5(10  fS 

P{S,V )  =  Pe(»0+- — -  4>(5)dS 

•  In 


F(’/)  r  . 

PAV)  +  — [E{S,V)  -  EAV )} 
V 


where  Eas.  (11)  and  (13)  have  been  used.  Eq.  (15)  is  the  usual  Mie-Gruneisen 
equation  of  state  and  does  not  require  knowledge  of  the  entropy  through  4>(5) . 


Equation  (15)  could  h-.ve  been  obtained  by  direct  integration  of  the  first 
form  of  Eq.  (5).  However,  this  method  shows  what  information  is  needed 
to  form  a  complete  equation  of  state  and  how  the  £-P«K  equation  does  not 
depend  upon  entropy. 

Consider  now  the  specific  heat  at  constant  volume  as  a  function  of 
T  and  V.  The  change  in  Cy(V,T)  along  an  isentrope  is  given  by 


/»- r\  .  fl£A  t/!£A/3E\ 

V  3V'/S  \  w/,  V  ztJv[wJ. 


using  Cv  =  T(B5/3r)j,  and  (Br/dE)5*  -(3P/BS)y.  Furthermore,  since 
T  =  (3£/3S)„  and  (3 P/?-T)v  =  (3S/3K)t, 


«). 


*  T 


(3P/3£), 

(3T/3E), 


*(*s\]  CCv)  {BP/B£'v 


3T\3K/t|  \3r/  (35/3E) 


.JM]  /^W3p\ 

_3rWvJ  \  oTj  vV>eJv 


(17) 


3  /bp\  1  ^cv/3r\ 

=  TCv[3fW/J  =  ~\dT/t 


(18) 


If  F  is  a  function  of  volume  only,  Eq.  (18)  shows  that  Gy(Y,T)  is  constant 
along  an  isentrope.  Thus  Cy(Y,T )  =  Cy[^p(S)]  is  some  function  of  the  entropy 
only.  Eq.  (9)  implies  that  when  T  is  a  function  of  volume  only. 


5  = 


(19) 
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Thus  for  F  «  T(K)  the  specific  heat  at  constant  volume  is  given  by 


(6( Y)\ 

Cv(V,T)  -g^— J 


One  special  case  of  interest  is  the  Debye  model  of  a  solid.  A  result 
of  this  model  is  that  the  specific  heat  at  constant  volume  is  a  function 
of  9d(V)/T  where  0D(V)  is  the  Debye  characteristic  temperature.  Identify¬ 
ing  0D(V)  from  the  Debye  model  with  0(V)  in  Eq.  (20)  and  using  Eq.  (10) 
yields  a  relationship  between  Gruneisen’s  ratio  P(V')  for  a  Debye  solid 
with  the  Debye  temperature  t9fl(’/) 


rm  = 


d  In  0D(V) 
d  In  V 


Slater  (Ref.  1>  suggested  a  relationship  between  r(K)  and  the  isen- 

tropes  of  the  system  based  on  the  Debye  model  and  the  assumption  that  in 

an  elastic  solid  the  ratio  of  the  longitudinal  sound  speed  to  the  trans¬ 
verse  or  shear  sound  speed  is  independent  of  volume.  This  assumption  is 

equivalent  to  assuming  that  the  ratio  of  the  isentropic  bulk  modulus  to 
the  isentropic  rigidity  modulus  is  volume  independent  or  that  Poisson's 
ratio  is  volume  independent.  The  Debye  characteristic  temperature  is 
given  by  (Ref.  1): 


0D{V)  =  A 


r  2V  I" 1/3 

Fs  «  c’J 


where  CL  and  CT  are  the  longitudinal  and  transverse  sound  wave  velocities 
respectively  and  A  is  a  constant. 

In  0D  =  In  A  In  2V  +  In  C£  +  “  In  (2  +  C\/C\ )  (23) 


d  In  Bq 
d  In  V 


j  d  In  CL 
3  +  d  In  V 


where  it  has  been  assumed  that  Cj/  CL  is  independent  of  volume, 
longitudinal  elastic  wave  velocity  is  given  by 

r  4 

CL  =  [V{K  +  -  #1)3 
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where  K  is  the  isentropic  bulk  modulus  and  fx  is  the  rigidity  modulus. 


UCL  -  i  I„  V  *  i  In  K  *  j  1„  (l  +  (26) 

Since  the  assumption  CT/CL  is  independent  of  volume  implies  fx/K  is  inde¬ 
pendent  of  volume, 


d  ln  CL  11  d  In  K 

d  In  V  "  2  2  d  In  V 

Now 


d  In  CL  (d2P/dV2  )s 

d  In  V  "  +  2  (dP/dV)s 

Combining  Eqs.  (29)  and  (24)  yields 


d  In  dD  v  (B2P/BV2)s  2 

r(V)  - - -  =  - - -  -- 

d  In  V  2  (BP/BV)  3 


(27) 


(28) 


(29) 


(30) 


This  is  the  Slater  relation  for  Gruneisen’s  ratio  in  terms  of  the  deriva¬ 
tives  of  the  isentropes.  It  has  been  derived  on  the  basis  of  the  Debye 
model  and  the  assumption  that  the  ratio  of  the  sound  speeds  is  independent 
of  volume.  Since  F  has  been  assumed  dependent  or.  volume  only,  the  ratio 
of  the  partial  derivatives  on  the  right  hand  side  of  Eq.  (30)  is  also 
dependent  on  volume  only,  that  is,  independent  of  entropy.  Since  the 
relation  must  hold  then  on  any  isentrope  at  a  given  volume,  it  must  hold 
on  the  isentrope  S  =  0  which  is  also  the  isotherm  T  =  0. 


Dugdale  and  MacDonald  (Ref.  2)  propose  that  the  Slater  relation  be 
modified  to 


V  B2(P’/2/3)/BF2  1 

2  B(PV'2/3)/Bl'  "1 


(31) 
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where  again  P  refers  to  the  isotherm  T  =  0 °K.  The  Dugda le-MacDonald 
relation  has  been  used  in  extensive  calculations  by  Rice,  McQueen,  and 
Walsh  (Ref.  3)  based  on  shock  data. 

In  summary,  the  usual  form  of  the  Mie-Gruneisen  equation  of  state 
is  that  of  Eq.  (15),  which  is  a  consequence  of  assuming  F  =  T ( l'* ) .  If  the 
Debye  model  of  a  solid  is  used,  F(P)  and  the  Debye  characteristic  tempera¬ 
ture  are  related  by  Eq.  (21).-  Substitution  of  the  expression  for  the 
Debye  temperature,  Eq.  (22),  into  Eq.  (21)  plus  the  assumption  of  constant 
Poisson’s  ratio,  yields  the  Slater  relation  between  F(  P).  and  the  0°K 
isotherm. 

Equation  (15)  may  be  written 

F(P) 

P  -  P  ( V )  =  -  [£  -  E  <P)]  (15) 

*  P 

The  point  P,  V,  E  refers  to  any  arbitrary  equilibrium  state  in  the 
region  of  the  pressure- volume  plane  for  which  Eq.  (15)  is  valid.  P^(P ) 
and  £^(P)  are  the  compressional  pressure  and  energy  at  that  same  volume. 
Equation  (15)  says,  then,  that  the  thermal  component  of  the  pressure  at 
volume  P  is  equal  to  F(P)/P,  a  function  of  volume,  times  the  thermal  com¬ 
ponent  of  the  energy  at  the  same  volume.  In  particular,  if  the  pressure- 
volume,  energy  state  under  consideration  lies  on  a  Hugoniot  curve  centered 
or  state  P0,  PQ,  £q,  Eq.  (15)  becomes 

PB  -  Pc(F)  =  [Eb  -  £c(P)]  (32) 

Subtraction  of  Eq.  (32)  from  Eq.  (15)  yields 

F(P)  r 

P-PB  £g]  (33) 

Equation  (33)  now  says  that  at  volume  P  the  difference  in  thermal  pressure 
between  the  state  P,  P.  £  and  the  Hugoniot  centered  at  PQ,  PQ,  £0  is  equal 
to  F(P)/P  times  the  difference  in  thermal  energy  between  state  P,  P.  £  and 
the  Hugoniot.  Since  0°K  isotherms  are  not  directly  measurable  as  arc 
Hugoniot  curves,  it  is  convenient  to  use  the  equation  of  state  in  the  form 
of  Eq.  (33)  with  the  Hugoniot  as  a  reference.  If  F(P)  is  assumed  to  have 
a  constant  value  G  and  the  experimental  Hugoniot  is  expressed  as  a  function 
of  p.  =  {p/p0)-  1,  Eq.  (33)  can  be  written 
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P 


PH{p)  +  GpE' 


(34) 


where  o  -  \/V  and  E' 
The  Hugoniot  internal 


=  E  ~  Eh,  the  thermal  energy  offset  from 
energy  and  pressure  are  related  through 


the  Hugcniot. 


-  £o+7W 

£°  +  2  p 


-  V) 


Substitution  of  Eq.  (35)  into  (34)  yields 


(35) 


P 


-y]  +  cp(£  -  £o> 


(36) 


where  E  is  the  total  internal  energy  and  £0  is  the  internal  energy  at  the 
foot  of  the  Hugoniot  Pg{p).  Equation  (36)  is  valid  for  G( ^  as  well  as 
G  constant-  However,  if  C  or  T  were  to  have  an  energy  or  temperature 
dependence,  the  P-V-E  equation  of  state  would  not  have  had  the  form  of 
Eq.  (15)  which  led  to  Eq.  (36).  If  f  is  a  function  of  the  energy  as  well 
as  volume,  the  generalized  Mie-Gruneisen  equation  of  state  is  given  by 
integration  of  the  differential  equation  (5)  defining  F(£,  V).  Taking 
E  and  V  as  independent  variables 


ns,  V ) 


(5) 


Integrating  at  constant  V 

P  ~  P  (V)  «  —  f  r{E.V),dE  (37) 

r  V 

where  £^(F)  and  Pf(V)  are  the  energy  and  pressure  at  some  reference  point 
at  the  volume  V. 

3.  SHOCK  WAVE  MEASUREMENTS  OF  EQUATION  OF  STATE 

A  series  of  experiments,  in  which  shock  waves  of  varying  intensity- 
are  produced  in  a  material  and  measurements  of  the  resulting  shock  and 
particle  velocities  are  made,  yields  a  -ingle  Hugoniot  curve  uniquely 
specified  by  the  initial  state.  The  specific  volume  V  and  pressure  P 
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are  related  to  the  shock  and  particle  velocities  through  the  equations 


V 


V 


0 


(38) 


and 


P  =  P0Uu 

which  express  the  conservation  of  mass  and  momentum,  respectively.  The 
energy  along  the  Hugoniot  curve  is  given  by 

E  =  Eo  *7  Wo  -  V)  (40) 

where  E0  is  the  energy  at  P  =  0,  V  =  V..  As  stated  earlier,  a  single 
Hugoniot  curve  is  not  an  equation  of  state:  it  is  only  a  particular  pressure- 
volume  curve  along  which  the  energy  can  be  computed  from  Eq.  (40).  A 
P-V-E  equation  of  state  consists  of  knowing  the  energy  over  a  region  of 
the  P~V  plane.  For  each  initial  state  there  exists  a  different  Hugoniot 
curve.  Therefore,  by  varying  the  initial  state  of  a  system,  shock-wave 

experiments  can  be  used  to  generate  many  Hugoniot  curves  in  the  P~Y  plane 
which  then  may  be  used  to  infer  the  dependence  of  the  internal  energy  on 
volume  and  pressure- 

In  the  work  to  be  described,  preheated  material,  porous  material,  and 
multiple  shocks  have  been  used  to  vary  the  initial  state.  The  preheating 
technique  consists  of  initially  heating  the  material,  aluminum  in  this 
case,  at  zero  pressure  to  near  its  melting  point.  The  initial  volume 
increases  cue  to  thermal  expansion  and  the  initial  energy  increases  by 
an  amount  A E  -  fCpdT.  The  double  shocking  technique  consists  of  reflecting 
the  initial  shock  induced  in  the  material  from  a  material  of  high  shock 
impedance,  tungsten  for  example.  The  state  behind  the  first  shock  then 
serves  as  an  initial  state  for  the  reflected  second  shock.  The  use  of 
porous  materials,  described  first  in  the  open  literature  by  the  Russians 
(Ref.  4),  produces  a  large  increase  of  thermal  internal  energy  in  the 
sample  during  shock  compression  due  to  the  irreversible  collapse  of  the 
porous  structure.  Due  to  the  large  thermal  energies,  Hugoniot  curves  of 
both  positive  and  negative  slopes  are  observed.  (Refs.  5,  6,  7)  Such 
behavior  of  the  Hugoniot  curves  of  metals  can  be  explained  by  using  the 
ordinary  Mie-Gruneisen  equation  of  state  with  i  =  1  (F)  or  a  more  elaborate 
equation  c.f  state  taking  account  of  the  thermal  contribution  of  the  conduction 
electrons. 
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SECTION  III 


EXPERIMENTAL  TECHNIQUE 


To  obtain  Hugoniot  points  for  aluminum  which  do  not  all  lie  on  a 
single  curve  the  three  approaches  described  in  Section  II  of  this  report 
have  been  attempted.  These  are  single  shocking  porous  material,  double 
shocking  porous  material,  and  preheating  the  material;  only  the  first 
method  was  used  extensively.  In  studying  Teflon,  the  first  method  was 
used  exclusively.  In  this  portion  of  the  report  are  described  sample 
preparation  and  experimental  technique. 

i .  METHOD  OF  OBTAINING  HUGONIOT  DATA 

Hugoniot  points  were  determined  by  the  impedance  match  method 
(Ref.  3),  using  2024  aluminum  as  a  standard,  as  illustrated  below  for  a 
porous  aluminum  specimen.  At  the  driver-porous  aluminum  interface  a  shock 
is  transmitted  into  the  porous  aluminum  and  a  rarefaction  is  reflected 
back  into  the  2024  driver.  Since 
the  Hugoniot  and  relief  isentrope 
cross  curves  in  the  pressure  (P)- 
particle  velocity  (u)  plane  for 
2024  aluminum  are  known, 

(Refs.  3,8)  a  measurement  of  the 
free-surface  velocity  of  the 
2024  driver  and  the  shock  veloc¬ 
ity  through  the  sample  suffice  to 
fix  a  point  on  the  Hugoniot  of 
the  sample.  The  measurement  of 
the  driver  free-surface  veiocitv 
gives  the  point  C  on  the  P-u  dia¬ 
gram  of  figure  1,  indicating  that 
the  initial  shock  state  in  the 
driver  is  state  B.  All  of  the 
states  which  the  2024  can  reach 
from  state  5  through  a  reflected 


FIG.  1  IMPEDANCE  MATCH  METHOD  OF 
DETERMINING  A  HUGONIOT 
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isentropic  rarefaction  lie  on  the  cross  surve  CDB.  Since  there  must  be 
continuity  of  pressure  and  particle  velocity  at  the  2024-porous  aluminum 
interface,  the  shocked  state  in  the  porous  aluminum  must  also  lie  on  this 
cross  curve.  Equation  (39)  requi res  that  the  Hugoniot  point  of  the  porous 
material  lies  on  a  ray  through  the  origin  of  slope  p^V/m,  where  m  is  the 
porosity,  p0  is  the  crystal  density  of  aluminum,  and  U  is  the  shock  veloc¬ 
ity  through  the  sample.  Hence  the  intersection  of  this  ray,  OD,  with  the 
cross  curve,  CDB,  yields  a  point  on  the  Hugoniot  of  porous  aluminum. 

2.  DESCRIPTION  OF  EXPERIMENTS 
a.  Rc om  Temperature  Shots 

Shocks  were  induced  into  2024  aluminum  drivers,  upon  which  the  samples 
were  mounted,  by  high  explosive  in  direct  contact  or  by  the  impact  of  an  ex¬ 
plosively  accelerated  flying  plate.  The  direct  contact  driver  system  con¬ 
sisted  of  an  8- inch-di ameter  plane  wave  lens2  ini tiating a  pad  of  9404  HE, 

4  inches  long  and  8  inches  in  diameter,  in  contact  with  a  3/8-  or  1/2-inch- 
thick  2024  aluminum  driver  plate. 

Most  of  the  experiments  required  flying-plate  systems  because  of  the  high 
pressures  desired.  The  samples  were  mounted  on  a  3/8-  or  1/4-inch-thick  2024 
aluminum  driver  plate.  A  strong  shock  wave  was  produced  in  the  driver  by 
impact  of  a  1/8-  or  3/16-inch- thick  stainless  steel  or  brass  flyer  plate 
(see  figure  2).  In  all  cases  an  8- to- 9- inch-di ameter  explosive  pad,  ini¬ 
tiated  by  a  plane  wave  lens,  was  used  to  accelerate  the  flyer  plate  across  a 
3/4-ineh  or  1-inch  air  gap.  To  prevent  spalling  and  breakup  of  the  flyer, 
it  was  separated  from  the  explosive  pad  by  a  1/16-inch  Plexiglas  buffer. 

All  velocities  were  determined  by  optical  techniques  with  a  rotating 
mirror  streak  camera  at  a  writing  speed  of  3.73  mm/jusec.  Multiple  slits 
were  used  and  illumination  was  by  an  explosive  argon  light  source.  Mea¬ 
surement  of  shock  velocity  through  the  samples  was  made  by  small  glass 
mirrors  bonded  to  the  driver  free  surface  end  aluminized  Mylar  bonded  to 
the  steel  rings  in  which  the  porous  samples  were  pressed.  Destruction  of 
the  mirrors  on  the  driver  surface  gave  the  arrival  time  of  the  shock  at 
the  driver-sample  interface  and  the  destruction  of  the  Mylar  gave  the 
arrival  time  of  the  shock  at  the  sample  free  surface. 

Measurement  of  the  free-surface  velocity  of  the  2024  aluminum  driver 
plate  is  made  by  mounting  a  mirror  inclined  at  a  small  angle  on  the  driver 

2  P-80,  P«ntf=. 
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FIG.  2  FLYING  PLATE  DRIVER  ASSEMBLY 

surface  (Ref.  9).  The  mirror  is  mounted  so  that  the  slit,  the  normal  to 
the  driver  surface,  and  the  normal  to  the  inclined  mirror  ail  lie  in  the 
same  plane,  as  shown  in  figure  3.  The  motion  of  the  free  surface  destroys 
the  inclined  mirror  continuously  over  a  small  time  interval  giving  an 
angular  cutoff  of  reflected  light  intensity  on  the  film  (see  figure  4b). 
The  mirror  angle  is  chosen  so  that  the  velocity  of  the  point  of  collision 
is  supersonic,  thereby  eliminating  the  possibility  of  a  premature  cutcff 
due  to  stress  waves  in  th«  mirror.  The  aluminized  surface  of  the  mirror 
is  covered  by  a  thin  (1  mil)  steel  shitn  to  prevent  premature  cutoff  by 

DIRECTION  OF  CAMERA  SLIT 


STEEL  SHIM. 


SHOCK  ARRIVAL 
MIRROR  ■ 


INCLINED  MIRROR 


SHOCK  ARRIVAL 
MIRROR  . 
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t  f  y  2024  -  ALUMINUM  DRIVER  PLATE*  Tf  >V 
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FIG.  3  INCLINED  MIRROR  TECHNIQUE  FOR  FREE- 
SURFACE  VELOCITY  MEASUREMENT 
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microjetting  or  by  a  shock  in  the  ambient  gas.  The  mirrors  mounted  on  the 
surface  on  either  side  of  the  inclined  mirror  indicate  shock  arrival  at 
the  free  surface  and  any  tilt  the  shock  front  may  have.  Knowledge  of  the 
initial  angle  at  which  the  inclined  mirror  is  mounted,  the  writing  speed 
of  the  camera,  the  magnification  of  the  camera  permit  calculation  of  the 
free-surface  velocity  from  the  angle  of  light  cutoff  on  the  film. 

Since  the  plane  area  on  the  driver  having  uniform  pressure  and  simul¬ 
taneous  breakout  time  was  less  than  4  inches  in  diameter,  much  care  was 
used  in  layout  of  samples  and  free-surface  mirrors.  Samples  were  mounted 
at  the  same  distance  from  shot  center  as  the  driver-state  inclined  mirrors 
in  order  to  minimize  the  effect  of  radial  pressure  variation.  A  typical 
streak  camera  record  is  shown  in  figure  4.  Figure  4a  is  a  still  picture 
of  the  shot  face  showing  the  mirrors  with  an  image  of  the  three  streak 
camera  slits  exposed  over  them.  Slits  I  and  3  record  data  from  the  two 
samples  being  studied.  On  slit  1  (slit  3  is  identical),  mirrors  1  and  2 
are  on  the  driver  and  record  the  arrival  time  of  the  incident  shock  at  the 
driver- sample  interface.  The  aluminized  Mylar  on  the  circular  sample  re¬ 
cord?  the  arrival  of  the  shock  at  the  sample  free  surface.  Mirror  6  is  an 
inclined  mirror  mounted  to  measure  the  sample  free-surface  velocity  in  the 
manner  just  described.  The  free-surface  velocity  of  the  aluminum  driver 
is  measured  on  slit  2  by  inclined  mirrors  ?  and  8  while  mirrors  3,  4,  and 
5  are  on  the  driver  surface  to  detect  shock  arrival  and  tilt. 

The  streak  camera  record  from  this  shot  is  shown  in  figure  4b.  From 
slit  1  (and  similarly  for  slit  3)  the  transit  time  At  of  the  shock  through 
the  sample  and  the  inclined  mirror  cutoff  yielding  sample  free-surface 
velocity  arc  obtained.  From  slit  2  the  inclined  mirror  cutoffs  yielding 
driver  free-surface  velocity  are  determined. 

To  reduce  extraneous  shock  light  a  methane  atmosphere  surrounded  the 
samples.  This  proved  particularly  effective  on  multiple  slit  shots. 

Similar  techniques  were  used  in  obtaining  high  pressure  data  by  double 
shocking  of  initially  porous  aluminum.  A  small  tungsten  block  was  placed 
over  half  of  the  sample  and  two  slits  observed  the  sample.  One  of  the 
slits  recorded  the  shock  velocity  through  the  sample  in  the  manner  just 
described  while  the  other  slit,  observing  the  half  of  the  sample  covered 
by  tungsten,  recorded  the  free-surface  velocity  of  the  tungsten.  A  third 
slit  recorded  the  2024  aluminum  driver  free-surface  velocity  just  as  did 
slit  2  in  figure  4.  The  free-surface  velocity  of  the  driver  and  the  shock 
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velocity  through  the  sample  yield  the  state  behind  the  initial  shock  in 
the  sample  just  as  before.  Upon  reaching  the  sample- tungsten  interface, 
the  initial  shock  is  transmitted  as  a  shock  into  the  tungsten  and  re¬ 
flected  as  a  second  shock  back  into  the  sample.  A  measurement  of  the 
free-surface  velocity  of  the  tungsten  together  with  the  Hugoniot  of  tung¬ 
sten  as  determined  by  McQueen  and  Marsh  (Ref.  10)  yield  the  state  behind 
the  transmitted  shock  in  the  tungsten.  Because  of  continuity  of  pressure 
and  particle  velocity  across  an  interface,  the  same  pressure-particle 
velocity  state  exists  behind  the  reflected  wave  in  the  sample.  The  final 
specific  volume  and  energy  behind  the  reflected  wave  may  then  be  calculated 
by  appropriately  applying  the  jump  conditions  across  the  reflected  shock. 

b.  Experiments  usihg  Preheated  1060  Aluminum 

Solid  1060  aluminum  samples  were  heated  to  temperature  in  excess  of 
500°C  and  then  impacted  by  an  explosively  accelerated  brass  plate.  Mea¬ 
surements  of  shock  velocity  through  the  sample  and  of  the  impact  velocity 
suffice  to  determine  a  Hugoniot  point  for  the  sample  if  the  Hugoniot  of 
the  flyer  plate  is  known.  The  experimental  arrangement  for  studying 
heated  samples  is  shown  in  figure  5-  Two  samples  of  1060  aluminum,  one 
to  be  heated  and  the  other  as  a  cold  control  sample,  are  mounted  in  an 
aluminum  base  plate  such  that  the  flyer  plate  strikes  the  samples  directly. 
The  sample  to  be  heated  (on  the  left  in  figure  5)  is  supported  in  a  quartz 
enclosure  to  insulate  it  thermally  from  the  base  plate.  It  is  heated  by 
a  resistance- wi re  heating  element  and  the  temperatures  of  both  the  hot  and 
cold  samples  are  monitored  by  thermocouples  as  shown.  The  face  of  the  shot 
is  viewed  through  the  slits  of  the  streak  camera.  Using  the  polished 
surfaces  of  the  samples  to  reflect  light  from  an  argon  light  source,  the 
shock  velocity  is  determined  by  recording  the  time  of  arrival  of  the  shock 
at  two  thicknesses  on  the  sample  in  much  the  same  manner  as  with  the  porous 
samples.  The  free-surface  velocity  is  determined  by  measuring  its  time  of 
flight  across  the  central  gap  to  a  quartz  mirror. 

Independent  measurements  of  the  flyer  plate  velocity  are  made  by  two 
sets  of  mirrors  mounted  in  the  flight  space  on  the  same  radius  as  the 
samples.  The  successive  destruction  of  the  reflectivity  of  the  central 
mirror  and  then  the  two  outer  mirrors  of  each  set  gives  the  time  of  flight 
over  a  5/32-ir.uh  gap  just  prior  to  flyer-sample  impact. 
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FIG.  5  EXPERIMENTAL  ASSEMBLY  FOR  FIRING  PREHEATED  SAMPLE 


Directly  behind  the  base  plate  is  a  0. 001-inch-thick  brass  or  stain¬ 
less  steel  radiation  shield.  It  was  found  desirable  in  two  of  the  shots 
include  this  heat  shield  to  prevent  excessive  heating  of  the  flyer  plate 
and  the  explosive  pad.  The  flyer  plate  of  1/8-inch  brass  was  accelerated 
over  a  1-inch  free  run.  Brass  was  used  in  these  experiments  since  exten¬ 
sive  Hugoniot  data  on  brass  are  available  from  the  work  of  McQueen  and 
Marsh  (Ref.  10).  Prior  to  firing,  the  entire  assembly  was  pumped  down  to 
pressures  below  5  £*Hg. 
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3.  SAMPLE  PREPARA i  ON 


To  prepare  the  samples,  aluminum  powder  of  high  purity  was  pressed 
in  steel  rings  to  the  desired  density.  The  values  of  porosi'v  a,  defined 
by  a  =  p0/p<a  where  pQ  is  crystal  density  and  p^  is  the  apparent  density, 
were  set  at  1.4,  1,7,  and  2.0. 

The  commercial  aluminum  powder  (Reynolds  ff  40  Atomized  Powder)  had 
an  average  particle  size  of  30  ±  5  microns  but  the  rather  lar^e  langv*-  of 
particle  sizes  made  sieving  necessary.  The  samples  of  porosity  a  1.4 
and  a  =  1.7  were  made  from  powder  which  was  passed  through  a  ff  65  mesh 
sieve  (208-micron  openings).  In  order  to  press  samples  of  porosity 
r  =  2.0  that  were  strong  enough  to  be  handled,  it  was  necessary  to  use  a 
larger  particle  size.  These  samples  were  made  from  powder  which  passed 
through  a  ff  50  mesh  sieve  (300-micron  openings)  but  not  through  a  ff  65 
mesh  sieve.  None  of  the  samples  required  sintering. 

No  experiments  have  been  performed  on  this  project  to  evaluate  the 
effect  of  particle  size  on  shock  wave  propagation  in  porous  aluminum. 
Russian  workers  (Ref.  6)  made  such  studies  on  porous  copper  and  lead  and 
concluded  that,  in  the  100-300  kbar  range,  changing  the  particle  size 
from  0.5  g  to  100  f*  did  not  affect  the  shock  velocity.3  It  would  seem 
that  particle  size  would  be  a  significant  variable  in  the  low  pressure 
region  where  compaction  of  the  porous  structure  is  incomplete.  However, 
no  experiments  have  been  performed  to  evaluate  the  effect  of  particle 
size  on  wave  shape  or  on  the  approach  to  equilibrium  in  the  pressure 
range  of  interest  here.  Both  theoretical  and  experimental  studies  of 
this  problem  are  needed  to  establish  the  validity  of  using  porous  materials 
to  attain  high  temperature- low  compression  states  by  shock  wave  techniques. 

To  make  a  sample,  a  weighed-portion  of  powder  was  carefully  leveled 
in  a  massive  steel  piston  cylinder  arrangement  which  held  the  steel  support 
ring.  After  pressing  several  times  to  a  thickness  stop,  the  sample  was  re¬ 
moved  in  its  ring  and  weighed  and  measured  for  density  and  thickness.  The 
uncertainty  in  the  average  density  of  a  typical  sample  was~0.  15  percent. 
Optical  records  of  shock  transit  times  through  the  specimens  indicated 
that  they  were  of  uniform  density.  A  thickness-to-diameter  ratio  for  the 
samples  was  chosen  to  provide  a  central  region  of  uniaxial  strain  (free 


However,  for  the  large  particles,  the  shock  rise  ti*e  say  be  no  longer  negligible,  and  thus  the  qaestion 
say  arise  as  to  shat  portion  of  the  profile  is  being  detected  as  "shock  arrival.** 
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Table  I 

TEFlJON  SAMPLE  measurements 


from  edge  effects)  sufficiently  large 
to  perform  precise  measurements.  A 
typical  sample  was  3 '16-inch  thick  and 
1-1/2  inches  in  diameter,  surrounded  by 
a  steel  ring  (1/8-inch  wall  thickness) 
and  covered  on  one  side  with  0.007- inch- 
thick  aluminized  Mylar. 

Solid  aluminum  samples  were  simply 
machined  and  lapped  from  1060  aluminum. 

The  solid  and  porous  Teflon  samples 
were  obtained  from  Avco  Corporation 
(for  description  of  preparation,  see 
Ref.  11),  in  the  form  of  disks  nominally 
2  inches  in  diameter  ai:d  1/2- inch  thick. 

They  were  sawed,  ground,  and  dry  lapped 
to  size  (nominally  3/16-inch  thick)  and 

their  porosity  was  determined  by  direct  measurement  of  dimensions  and 
weights.  Teflon  sample  densities  are  shown  in  Table  I. 


SAMPLES0 

THICKNESS 
( i aches  ) 

INITIAL  DENSITY 
(g/cn3) 

2-A 

0. 1870 

2. 1639  t  0  0046r 

2-B 

0. 1859 

2.1697  t  0.0035' 

6 

0. 4949? 

1.472 

10 

0. 4944 

1.489 

10-A 

0.1874 

1 . 525 

10-B 

0. 1909. 

1.513 

3 

~0. 49117 

1.511 

9 

~0  4926 

1.437 

1-A 

0. 1885 

0.733 

1-B 

0. 1742 

0.778 

а.  Slices  of  as -received  staples  indicated 
by  letter  on  sasrple  nonber 

б.  ftisetisiocs  of  escb  as-received  sasple 
•tried  by  several  nils.  Average 
values  are  given. 

c.  A  valne  of  -  2.1668  g/cn?  «as  as--d 
in  calculating  the  porosity  a. 
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SECT! OS  IV 


EXPERIMENTAL  RESULTS 


1.  ALUMINUM  DATA 

The  measured  free-surface  velocities  cf  the  2024  aluminum  driver  plates 
and  the  measured  shock  and  free-surface  velocities  of  the  solid  and  porous 
aluminum  samples  are  presented  in  Table  II.  The  pressure-particle  velocity 


Table  II 

SOLID  AND  POBQUS  ALUMINUM  SHOCK  DATA 


2024  DRIVER 

SPECIMEN'9 

SHOT  SO- 

FREE-SURFACE 

VELOCITY 

(aa/fixec) 

Iaitial 
IX.‘Bsa  t  jr 
l*/c«3) 

Shock 

Velocity 

(aa/pxec) 

Free-Sarfwce 

Velocity 

(n/piec) 

11.470 

2.89 

a  *  1.0 

2.70 

7.31 

2.97 

11.570 

3.6C 

2.70 

7.75 

— 

11.563 

5.03 

2.70 

8.77 

... 

11.460 

5.35 

2.70 

9.06 

5.42 

11.562 

5.61 

2.70 

9.13 

— 

11.477 

2.89 

a  =  1-4 

1.93 

5.03 

3.23* 

10.591 

5-07« 

1.92 

7-15 

5-43 

1C. 592 

5.86 

1.92 

7.79 

5.93 

11.155 

3-56 

=  =  1.7 

1.59 

5.43 

3.96* 

11.286 

3-96 

1.59 

5.82 

4.14* 

11.052 

4.02 

1.59 

5.93 

4.53 

10.376-1 

4.67 

1.59 

6.58 

5.21 

10.S76-2* 

4-67 

1.59 

6.50 

5.25 

10.926 

4.97 

1.53 

6.83 

5.53 

11.330 

5.04 

1.59 

6. SO 

5.52 

10.591 

5.07c 

1.5S 

6.91 

5-25* 

1G.894 

5.43 

1.59 

7.29 

5.95 

11.305 

5.72 

1.53 

7. 49 

-- 

10.592 

5.86 

1.58 

7.55 

6.24 

11.155 

3.56 

a  =  2.0 

1-35 

5.09 

4.51* 

11.286 

3.96 

1.35 

5.57 

4.16* 

11.062 

4.02 

1.35 

5.64 

5.11 

10.926 

4.97 

1.35 

6.72 

6.03 

11.330 

5.04 

1.35 

6.79 

5.84* 

10.894 

5.43 

i.35 

7.l6 

6.85 

11,379 

o«34 

1.35 

7.29 

-- 

11.305 

5.72 

1.35 

7.32 

— 

"All  were  itcuilh  3/26-  itch  thick,  except  fcr  the  3/t-iacS 

spcciaea  ei  Sk-t  Nc-  10.876-2. 

'  Alctim  shia  cx  free  ssrface- 
C  Adjwxted  ••!«  («ee  text}. 
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states  in  the  samples  vere  obtained  through  an  impedance  natch  solution 
using  the  2024  aluminum  as  a  standard,  as  described  in  Section  III.  The 
2024  aluminum  Hugoniot  and  release  cross  curves  were  obtained  from 
Dr.  !i.  G.  McQueen  of  the  Los  Alamos  Scientific  Laboratory.  The  processing 
of  the  porous  aluminum  data  consisted  of  first  plotting  the  specimen  shock 
velocity  against  the  driver  plate  pressure  that  was  inferred  from  the  mea¬ 
sured  free-surface  velocity  (using  the  LASL  cross  curves).  This  plot  of 
primary  data  (net  included  in  this  report)  served  as  an  initial  check  on 
the  scatter  in  the  data.  The  data  for  each  of  the  porosities  1.7  and  2.0 
defined  a  smooth  curve.  The  data  for  porosity  1.4  were  not  as  abundant  as 
for  the  other  porosities.  In  one  shot  (No.  10,591),  the  record  yielding 
the  free-surface  velocity  of  the  driver  plate  was  ambiguous.  Therefore, 
the  point  for  the  sample  of  porosity  1.7  was  made  to  lie  on  the  shock 
velocity-driver  pressure  curve  which  was  well  defined  by  the  other  data 
points.  In  so  doing,  the  point  for  the  1.4  porosity  sample  was  brought 
into  a  position  consistent  with  a  smooth  curve  through  the  remaining  points. 

The  pressures,  particle  velocities,  volumes,  and  energies  computed  from 
the  primary  data  by  use  of  the  impedance  match  technique  and  the  Rankine- 
Hugociot  jump  conditions  are  presented  in  Table  III.  Some  of  the  data 
reported  in  Tables  II  and  III  differ  slightly  from  those  reported  earlier 
in  the  progress  reports.  The  data  have  been  rechecked  for  errors  in  calcu¬ 
lations,  and  questionable  shot  records  have  been  reanalyzed.  The  present 
results  reflect  the  changes  arising  in  the  reexamination. 

The  porous  aluminum  data  are  plotted  on  a  shock  velocity-particle 
velocity  diagram  in  figure  6.  Plotted  on  the  same  figure  are  some  porous 
aluminum  data  obtained  by  Avco  (Ref.  11).  Straight  lines  represent  all  the 
data  quite  well  for  a  =2.0  and  a  =  1.7  over  the  rangt  investigated.  The 
three  points  for  a  =  1.4  lie  on  a  straight  line.  The  three  highest  pressure 
Avco  points  for  a  =  1.5  lie  between  the  lines  for  x  =  1.7  and  a  =  1.4  but 
indicate  a  different  slope.  The  two  lower  Avco  a  =  1.5  points  appear  ques¬ 
tionable  (indicated  by  parentheses)  as  they  imply  pressure-volume  states 
which  are  quite  scattered  from  the  remainder  of  the  data. 

Representation  of  shock  data  for  porous  material  by  a  linear  shock 
velocity-particle  velocity  relationship  is  invalid  for  the  purpose  of  ex¬ 
trapolation  to  lower  pressures.  The  relationship 

V  =  C  +  Su  (41) 
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Table  III 

SOLID  AND  POflOUS  ALUMINUM  HUGONIOT  DATA 


SHOT  NO. 

INITIAL 

VOLUME 

(cnVj) 

SHOCK 

VELOCITY 

(n/fut  c  ) 

PARTICLE 

VELOCITY 

(ss/^scc) 

PRESSURE 

(kbar) 

FINAL 

VOLUME 

{caVg) 

COMPRESSION3 
u  =  p/p0  -  1 

INTERNAL 
ENERGY 
CHANGE 
(10®  ergs/*} 

11,477 

|  | 

7,31 

1.45 

286 

0.297 

0.250 

10.65 

11,570 

0.37u 

7.75 

1.80 

377 

0.284 

0.302 

16.21 

11.567 

8.77 

2.46 

583 

0. 266 

0.391 

30.31 

11,460 

0.370 

9.06 

2.58 

633 

0.265 

0.3% 

33.23 

11,562 

0.370 

9.13 

2.71 

667 

0.260 

0.423 

36.56 

11,477 

0.519 

5.08 

1.88 

184 

0.328 

0.128 

17-57 

10,591 

0.520 

7.15 

3.03 

418 

0.299 

C.237 

46.18 

10,592 

0.521 

7.79 

3.42 

511 

0.293 

0.262 

58.25 

11.155 

0.629 

5.43 

2.43 

210 

0.347 

0.066 

29.61 

11.286 

0.631 

5.82 

2.67 

246 

0.341 

0.985 

35.67 

11.052 

0.629 

5.93 

2.69 

254 

0.344 

0.075 

36.19 

10,876-1 

0.630 

6.58 

3.05 

319 

0.338 

0.094 

46.57 

10.876-2 

0.629 

6.50 

3.06 

316 

0.333 

0.111 

46.76 

10.926 

0.631 

6.83 

3.20 

347 

0.335 

Q.1C4 

51.35 

11.330 

0.631 

6.90 

3.25 

355 

0.334 

0.107 

52.72 

10,591 

0.633 

6.91 

3.26 

356 

0.334 

0.107 

53.22 

10.894 

0.630 

7.29 

3.44 

397 

0.333 

0.111 

58.95 

11.305 

0.631 

7.49 

3.58 

426 

0.330 

0.121 

64.11 

10.592 

0.635 

7.55 

3.66 

436 

0.327 

0.131 

67.14 

11.155 

0.741 

5.09 

2.5B 

177 

0.365 

0.0136 

33.19 

11.286 

0.741 

5.57 

2.82 

212 

0.364 

0.0164 

33.85 

11.062 

0.741 

5-34 

2.85 

218 

0.367 

0.CG31 

40.76 

10,926 

0.741 

6.72 

3.37 

305 

0.369 

0.0027 

56.73 

11.330 

0.740 

6.79 

3.41 

313 

0.369 

0.0027 

58.06 

10.894 

0.741 

7.16 

3.60 

318 

0.362 

0.0054 

64.90 

11.379 

0-744 

7.29 

3.85 

377 

0.351 

0.0541 

74.08 

11.305 

0.741 

7.32 

3.78 

373 

0.360 

0.0277 

71.05 

Valees  of  fx  based  oa  crystal  dcasity  P0  ~  2.70  g/ea  . 
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where  C  and  5  are  constants,  implies  a  pressure - vo 1 ume  Hugoniot  given  by 


^0  ,r  . 

-  C2 [1  -  V/mV  J 

P  =  -- - (42) 

[1  "5(1  -  V/mVQ  )]2 

where  V0  and  pQ  are  the  initial  specific  volume  and  density,  respectively, 
of  the  solid  material.  The  pressure  vanishes  at  in  Eq.  (42),  rather 

than  at  Kg  as  it  must  in  a  model  that  assumes  zero  crushing  strength  of  a 
porous  soLid.  Although  the  crushing  strength  is  not  actually  zero,  the 
low  pressure  extrapolation  of  Eq.  (39)  indicates  a  much  higher  pressure 
and  greater  curvature  of  the  Hugoniot  than  is  consistent  with  very  low 
pressure  experiments  carried  out  in  this  laboratory  on  another  project 
(Hef.  12). 

The  data  have  been  plotted  in  the  pressure- volume  plane  in  figure  7- 
Shown  also  on  the  figure  are  the  Hugoniot  of  solid  aluminum  from  Al'tshuler 
and  the  0°K  isotherm  from  the  McCloskey  program  (Ref.  13).  Measurements  of 
the  density  of  the  solid  1G60  aluminum  yie  lded  a  value  of  2.70  g/cm3.  The 
Russians  quote  a  value  of  2.71  g/cm3  for  the  density  of  aluminum  used  in 
their  experiments.  To  make  their  Hugoniot  comparable  to  the  data  reported 
here  each  of  their  volumes  at  a  given  pressure  was  changed  by  a  factor 
2.71/2.70.  This  maintains  the  shape  of  the  curve  but  translates  it.  slightly 
toward  the  right,  i.e.,  larger  volumes.  During  the  course  of  studying 
aluminum,  it  was  also  found  that  the  Los  Alamos  Hugoniot  curve  for 
2024  aluminum  could  be  made  to  agree  quite  closely  with  the  Soviet  Hugoniot 
for  pure  aluminum  by  translating  the  2024  curve  to  the  right  by  a  factor  of 
2.78/2.71,  the  ratio  of  the  initial  densities.  Smooth  curves  have  been 
drawn  through  the  porous  aluminum  data.  These  curves  do  not  represent  ana¬ 
lytical  fits  to  the  data  nor  are  they  obtained  from  the  straight  lines  of 
figure  6.  They  have  been  extrapolated  to  zero  pressure  and  crystal  volume 
(dashed  portion),  so  as  to  be  consistent  with  a  negligible  crushing  strength. 
In  the  absence  of  data  at  very  low  pressure,  there  appears  to  be  no  a  priori 
justification  to  make  any  other  type  of  extrapolation. 

The  two  highest  pressure  points  on  the  Hugoniot  for  porosity  m  =  2  have 
not  been  weighted  in  drawing  the  curve  through  those  data.  The  evidence 
suggests  that  the  abrupt  change  in  slope  indicated  by  those  points  is  not  a 
real  property  of  the  aluminum  but  rather  an  effect  of  shock  attenuation. 
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The  attenuation  is  due  to  the  overtaking  of  the  shock  by  the  rarefaction 
from  the  back  surface  of  the  flyer  plate.  There  is  no  conclusive  proof 
that  this  is  the  case,  but  calculations  based  on  sound  speeds  in  solid 
aluminum  indicate  that,  in  the  high-pressure  porous  aluminum  experiments, 
the  margin  of  safety  for  obtaining  data  free  from  attenuation  effects  is 
quite  small.  In  figure  7  two  points  on  the  Hugoniot  of  porosity  n  =  1.7 
are  labeled  t  -  3/16  inches  and  t  =  3/8  inches,  where  the  t  refers  to  the 
sample  thickness.  Both  samples  were  mounted  on  the  same  type  of  driver 
plate.  The  driver  under  the  3/8-inch  sample  was  struck  by  a  3/16-inch 
flyer  plate  while  that  of  the  3/16  sample  was  struck  by  a  1/3-inch  flyer. 
Attenuation  will  be  most  likely  for  the  smallest  ratio  of  specimen  thick¬ 
ness  to  flyer  thickness,  i.e.,  for  the  thickness  sample  in  the  present 
case.  Indeed,  it  appears  that  attenuation  did  affect  the  thicker  sample. 
The  few  experiments  performed  on  this  project  to  attempt  to  evaluate  shock 
attenuation  effects  indicate  a  need  for  more  studies  in  shock  wave  propa¬ 
gation  in  porous  materials. 

The  uncertainties  in  each  of  the  points  in  figure  7  were  calculated 
by  a  computer  program  which  takes  into  account  uncertainties  in  measured 
wave  and  free-surface  velocities  and  in  the  initial  density  of  the  sample. 

According  to  Eq.  (15),  if  the  Gruneisen  ratio  is  dependent  only  upon 
volume,  it  represents  the  ratio  of  the  thermal  pressure  to  the  thermal 
energy  density 


rm 


P  -  Pe(V) 

v - 

E  ~  Ee(V) 


(15) 


where  the  subscript  “ th ”  means  thermal.  In  cases  where  Gruneisen’s  ratio 
is  a  function  of  volume  plus  another  independent  variable,  Eq.  (15)  is  not 
the  Mie-Gruneisen  equation  of  state.  However,  a  parameter  X  analogous  to 
the  Gruneisen  ratio  may  be  defined  by 


MEM 


P  -  Pr(V ) 

V/  - 

E  ~  Er{V) 


(43) 


where  P  {V)  is  some  reference  curve  in  the  P  -  V  plane  and  Er(V)  is  the 
energy  along  that  curve.  If  P  (V)  and  £  ( K)  are  the  0°K  pressure  and 
energy,  then  X  is  the  ratio  of  the  thermal  pressure  to  thermal  energy- 
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density.  In  the  limit  as  the  curve  P  ( V )  is  aDoroached  at.  mnstani;  volume, 
k(E,V)  approaches  P(E, K) .  since  Eq.  (41)  approaches  the  partial  derivative 
which  defines  T(£, F). 

In  earlier  quarterly  reports  (Ref.  14)  some  values  of  X  were  computed 

from  the  experimental  data  using  the  Soviet  Kugoniot  (Ref.  6)  for  solid 

aluminum  as  a  reference.  Plots  of  the  resulting  values  against  volume 
indicated  a  decrease  of  X  with  decreasing  volume  and  increasing  energy. 

It  was  pointed  out  that  the  uncertainties  in  X  could  be  as  large  as  20  per¬ 
cent.  It  now  appears  that  this  estimate  of  the  uncertainty  in  X  was  quite 
optimistic . 

As  mentioned  earlier,  the  data  have  recently  been  reanalyzed.  The 
positions  of  the  aluminum  Hugoniots  for  porosities  a  =  1.7  and  m  =  2.0 
have  not  been  significantly  altered.  However,  the  Hugoniot  for  a  =  1.4 
has  been  shifted  toward  the  solid  Hugoniot,  particularly  in  the  low  pres¬ 
sure  region,  due  to  the  addition  of  a  new  point  (Shot  11,477)  and  the 
rejection  of  Shot  10,590  (because  of  the  poor  quality  of  the  film  record). 

The  resulting  curves  are  those  of  figure  7.  The  values  of  X  at  a  given 

volume  were  substantially  reduced.  Since  individual  points  are  subject  to 
some  scatter  and  the  resulting  values  of  X  are  so  sensitive  to  the  movement 
of  a  Hugoniot  point  in  the  P  -  V  plane,  it  appears  reasonable  to  draw  a 
smooth  curve  through  the  data  points  in  the  P  **  V  plane  and  then  to  compute 
values  of  X  from  the  curve  rather  than  from  individual  points.  It  also 
appears  to  give  X  a  more  unique  value  if  the  0°K  isotherm  is  used  as  a 
reference  curve  instead  of  the  solid  Hugoniot.  In  this  case,  X  =  KPth/£tk- 
The  0°K  isotherm  from  the  McCloskey  program  was  used. 

A  plot  of  the  points  against  taken  from  the  Hugoniot  curves 

and  0°K  isotherm  of  figure  7  is  shown  in  figure  8-  The  curve  obtained 
from  Hugoniot  r  =  2  appears  to  have  positive  curvnture  while  the  curves 
for  «  ~  1.0,  1.4  and  1.7  appear  to  have  negative  curvature.  The  points 
for  m  -  1  cluster  near  the  origin  as  thermal  pressure  and  energy  are  com¬ 
paratively  low  along  this  Hugoniot. 

The  values  of  X  calculated  from  each  of  the  Hugonict  curves  for 
a  =  1.4,  1.7,  and  2.0  are  plotted  against  volume  in  figure  9  and  against 
thermal  energy  in  figure  9(b).  The  curves  for  the  Hugoniot  m  =  1  are 
values  computed  by  the  Russians  from  the  0°K  isotherm,  using  the  Slater 
relation  for  P.  It  must  be  emphasized  that  these  values  of  X  are  calcu¬ 
lated  from  the  curves  drawn  through  the  data  and  that  the  drawing  of 
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FIG.  9  EFFECTIVE  GRUNEISEN  RATIO  A  FOR  ALUMINUM  AT  CONSTANT  FOROSITY 


slightlv  different  curves  would  yield  quite  different  \  values.  Also,  tne 
curves  Sve  been  forced  to  pass  through  \  -  2.1  at  1  =  V  e  and  P  -  0.  ’  The 

Hugoniot.  for  a  -  1.7  is  the  most  well  defined  of  the  three  as  it  is  based 

on  no."1  data.  The  Hugoniot  a  =  1.4  is  the  most  questionable  as  it  is  based 
upon  only  three  data  points.  An  analysis  of  the  above  results  was  made  to 
see  if  the  dependence  on  volume  aid  on  thermal  energy  could  be  resolved  in 
a  simple  way.  To  achieve  this,  curves  of  A  as  a  function  of  volume  at  con¬ 
stant  thermal  energy  and  as  a  function  of  thermal  energy  at  constant  volume 
were  plotted.  The  results  are  shown  in  figures  10  and  11.  A  maximum  of 

three  points  was  available  for  eac1.  cure.  The  minima  in  the  A  i*s.  £tJj 

curves  are  due  to  the  low  values  of  A  obtained  from  the  a  =  1.4  Hugoniot. 

As  this  Hugoniot  is  not  well  defined,  no  significance  can  be  attached  to 
the  minima.  It  is  felt  that  carrying  this  analysis  further  is  not  war¬ 
ranted  with  the  present  quantity  of  data.  Obtaining  a  large  amount  of  data 
so  that  a  statistical  analysis  could  be  applied  appears  to  be  the  most 
meaningful  method  of  truly  studying  the  behavior  of  Gruneiser.'s  ratio  Iron 
shock-wave  data. 

Three  high  pressure  points  were  obtained  by  the  multiple  shock  tech¬ 
nique  using  tungsten  as  a  shock  reflector.  The  data  are  presented  in 
Table  IV.  The  values  of  A  are  calculated  from  the  final  shock  state.  The 
internal  energies  of  the  two  highest  pressure  points  are  in  excess  of  the 
sublimation  energy.  In  considering  the  pressure  and  energy  range  covered 
by  these  data  compared  with  the  single  shock  data  on  porous  aluminum,  the 
most  significant  result  is  the  fact  that  although  it  does  not  appear  pos¬ 
sible  to  definitely  establish  the  energy  and  volume  dependence  of  A,  the 
values  obtained  do  not  vary  widely  and  are  quite  close  to  the  zero  pressure 
value  obtained  from  thermodynamic  data. 

The  Hugoniot  data  obtained  from  the  experiments  in  which  preheated 
solid  1060  aluminum  samples  wei e  used  arc  presented  in  Table  V.  For  the 
purpose  of  studying  Gruneisen’s  ratio  as  calculated  from  thermal  pressure 
and  energy  offsets  from  the  0°K  isotherm,  these  data  are  not  useful.  The 
offsets  in  the  preheated  Hugoniot  points  from  the  room  temperature  Hugoniot 
of  solid  1060  aluminum  were  too  small  to  be  considered  significant. 


Ai  (P  =  0,  V  =  fj),  A  =  I  =  V{CP/(£)y  =  VBjS/Cy  =  2.1  fer  alsnicsa,  there  Bj  is  the  issthersal  bilk 
sodslss  and  3  is  the  eolsse  coefficieit  of  theraal  expacsioa. 
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FIG.  10  EFFECTIVE  GRUNEISEN  RATIO  A  vs.  SPECiFiC  VOLUME 
AT  CONSTANT  THERMAL  ENERGY 


FIG.  11  EFFECTIVE  GRUNEISEN  RATIO  A  «.  THERMAL  ENERGY 
AT  CONSTANT  VOLUME 


Table  IV 

DOUBLE  SHOCK  POROUS  ALUMINUM  DATA 


SHOT  XO. 

IXITlAL 

VOLLVE 

(caVf  ) 

FIRST  SHOCK 

SECOND  SHOCK 

TOTAL 

EXEP.GY  3 
CIO9  er^»/g) 

* 

Pr'isore 

(Ibirl 

Fia<  1 

Vo lose 

<c=3/g} 

EstrgT 

CLaage 

HO9  «:p/g) 

Presssre 
debar  ) 

Fia*  1 
Vo lose 

(caVg) 

Eaergy 

Ch» age 

CIO9  ergs/g) 

10.876 

0.630 

319 

0.333 

46.57 

934 

44.30 

92.53 

1.90 

11.305 

0.631 

426 

0.330 

64.11 

1.261 

0.261 

58.20 

123.97 

1.87 

11.305 

0.741 

373 

0.360 

71.05 

1.244 

0.282 

63.05 

135.77 

2.03 

c  £  =  ifj  +  +  3.66  *  109  erg*/g. 
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Table  V 

HL’OQNTCT  DATA  FOR  PREHEATED  AND  ROOM  TEMPERATURE 
SOLID  1060  ALUMINUM 


To 

(°CI 

SHOCK 

VELOCITY 

(ss/fi>e  c  ) 

PARTICLE 
VELOCITY 
(sa/fiscc  ) 

FREE-S13FACE 
VELOCITY 
(ss/^sec  ) 

P1 

(kbar ) 

rl 

(e=3/g) 

542 

3.82 

2.75 

5.47 

643 

0.264 

44.6  *  109 

17 

3.76 

2.77 

5.48 

656 

0.252 

38.7  *  109 

555 

8.56 

2.76 

5.27 

611 

0.262 

43.3  x  109 

17 

8.82 

2.70 

5.14 

647 

0.257 

36.6  *  109 

543 

8.24 

2.23 

4.66 

475 

0.282 

30.1  *  109 

17 

8.14 

2.21 

4.56 

486 

24.3  r  109 

"  Energy  isrreue  relative  to  state  P  -  0.  T—  17°C- 

As  a  consistency  check  some  Hugoniot  points  for  solid  room  temperature 
1060  aluminum  of  measured  initial  density  2.70  g/cn3  were  obtained  by  the 
same  technique  used  on  the  porous  samples.  The  points,  shown  on  figure  7, 
are  in  good  agreement  with  the  Soviet  Hugoniot3  for  solid  aluminum  at  low 
pressures — and  oniy  fair  at  the  highest  pressures. 

2.  TEFLON  DATA 

Ten  data  points  on  solid  and  porous  Teflon  (see  Table  VI)  were  ob¬ 
tained  using  the  same  explosive  driver  systems  and  impedance  match  tech¬ 
nique  (2024  alumiium  standard)  as  was  used  in  obtaining  the  porous  aluminum 


Tabic  VI 

SOLID  AND  POROUS  TEFLON  DATA 


SHOT  SO. 

'*0 

<g.c»3) 

m 

■ 

c 

iu/jtsec) 

P 

(kbar ) 

■ 

(sa/ptc) 

V 

(e*3/g  ) 

££ 

{109  ergs/g) 

11.570 

0.7828 

4.374 

103 

2.985 

0.4057 

44.87 

11.563 

9.7782 

2_ 

7# 

5.646 

175.5 

3.985 

0.3780 

79.58 

11,562 

0.7828 

6.086 

207.0 

4.335 

0.3676 

94.12 

11.570 

1.517 

4.909 

190.3 

2.555 

0.3160 

32.64 

11.563 

1.524 

i. 

42 

6.207 

319.0 

3.373 

0.2995 

56.87 

11.552 

1.517 

6.683 

371.0 

3.66 

0.2981 

66.97 

11.570 

2.169 

5.589 

269.0 

2.215 

0.2782 

24.58 

11.363 

i. 

00 

6.830 

435.0 

2.945 

43.37 

11.450 

7.105 

476.0 

3.080 

47.50 

11.562 

7.262 

504.0 

3.20 

51.35 

3  Tie  shift  ia  tke  Soviet  Begsa jot  described  earlier  i»  sell  vitlis  experiaextel  error. 


38 


data.  In  figure  12  the  shock  velocity  is  plotted  against  the  particle 
velocity.  The  data  obtained  by  the  Avco  Corporation  are  also  included 
on  this  graph  for  comparison.  Figure  13  is  a  plot  of  th.  Hugoniot  pres¬ 
sure  as  a  function  of  volume.  As  with  the  porous  aluminum  data,  smooth 
curves  have  been  drawn  through  the  points. 

Since  there  is  no  theory  at  the  present  time  as  to  the  equation  of 
state  of  a  plastic,  and  since  the  present  data  are  preliminary,  it  is 
felt  that  further  analysis  should  wait  until  more  data  are  obtained.  How¬ 
ever,  some  estimates  of  the  Gruneisen  ratio  can  be  made  on  the  basis  of 
the  data.  In  the  absence  of  a  0°K  isotherm  for  Teflon  a  few  values  of  \ 
were  ccmputed  from  the  curves  drawn  through  the  Hugoniots  of  porous  Teflon 
using  the  curve  drawn  through  the  Hugoniot  of  solid  Teflon  as  a  reference. 
The  results  of  these  calculations  are  presented  in  Table  VII.  An  estimate 
of  F0,  Gruneisen's  ratio  at  zero  pressure,  was  made  for  Teflon  based  upon 
available  zero  pressure  data.  From  the  thermodynamic  definition  of 
Gruneisen’s  ratio 

VBrB 

r  =  — —  (44) 

Cy 

where  BT  is  the  isothermal  bulk  modulus  and  is  the  thermal  expansion 
coefficient.  Values  of  /5  and  Cy  were  found  to  be  (Ref.  15) 

£  =  16.5  x  10‘$/WC 

Cv  =  0.25  cal/g  (at  73°C) 

Br  was  estimated  from  the  data  of  Bridgman  to  be  43  kbar  (Ref.  16).  In¬ 
serting  these  values  into  Eq.  (44)  yields  a  value  1  e  =  0.313. 
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FIG.  12  PARTICLE  VELOCITY  vs.  SHOCK  VELOCITY  FOR  TEFLON 
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FIG.  13  HUGONIOTS  FOR  TEFLON 


Table  VII 

CALCULATED  VALUES  OF  A  FOR  TEFLON 


0  Hufosiot  of  Teflon  a  -  1.0  used  at  a  reference  curve. 
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